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We present analytic analyses of the joint density of states (JDOS) and the elastic scattering 
susceptibility (ESS) for d-wave superconductors under the linear dispersion approximation. The 
JDOS and ESS diverge on the same curves which are found to be the envelope functions related 
to the contour of constant energy of the Bogoliubov quasiparticles. The relative amplitudes of 
the divergence on the envelope functions are derived. It is also found that certain octet vectors 
come close to the vertices of the envelope curves when the Dirac-cone anisotropy is strong, which is 
relevant to the d-wave cuprates. 
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The scanning tunneling spectroscopy(STS) has long 
been applied to exploration of the local single-particle 
spectrum of high- Tc cuprates — . Recently, modulations 
of the local density of states (LDOS) have been stud- 
ied by the Fourier transformed (FT) STS and analyzed 
within the quasiparticle interference (QPI) picture that 
attributes the QPI pattern to the elastic Bogoliubov 
quasiparticle scattering from impurities — i^. The char- 
acteristic wave vectors seen in the FT-LDOS are inter- 
preted by an effective octet modelS-^'' that proposes seven 
wave vectors connecting the tips of the banana-shaped 
contours of constant quasiparticle energy (CCE). Numer- 
ical calculations using the T-matrix formalism^^^ on the 
scattering from a weak single impurity were in agreement 
with the octet model. 

As is known, the angle resolved photoemission spec- 
troscopy (ARPES) and the STS are complementary to 
each other, both measuring the single-particle electronic 
structure with the former being in the momentum space 
and the latter in the real one. A more intimate relation 
between them has been noted recently by comparing the 
two-particle information extracted from these two mea- 
surementgS, namely the susceptibility of the supercon- 
ductor to impurity embodied by the FT-LDOS intensity 
and the joint density of states by evaluating the autocor- 
relation of the spectral function observed by ARPES. A 
remarkable matching between the peak positions in the 
autocorrelated ARPES data and the octet vectors im- 
plies likely that the JDOS has an intrinsic relationship 
with the QPI patterns. However, the inherent physi- 
cal reason of the above-mentioned matching is unclear. 
Moreover, to the best of our knowledge, there still lacks a 
simple theoretical relation between the FT-LDOS, which 
is closely related to the imaginary part of the autocor- 
relation of the Nanibu-Gor'kov Green's function in the 
weak scattering limit, and the JDOS. 

In this paper, we attempt to uncover the inherent 
physics of the above-mentioned intimate relation between 
the JDOS and the QPI pattern. In particular, an enve- 
lope scenario is proposed to analyze the divergent be- 



havior of the JDOS, which reveals its highly structured 
pattern. The matching of the octet wave vectors and the 
JDOS peak positions is also investigated. 

Let us first derive the basic formulas for both the FT- 
STS and the JDOS as well as their intrinsic relationship. 
The ESS for the weak nonmagnetic impurity associated 
with the FT-STS is defined as 

5(p,w) - --Im V[G(k + p,c^+)a3G(k,a.+ )]n, (1) 
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where 



,+ = 



uj+iO'^ and G is the Nambu-Gor'kov Green's 



function of unperturbed d-wave superconductors, 
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with ^k = tss^^^ — /i the normal-state dispersion 
measured with respect to the chemical potential, Ak = 
— 2Ao(cos kx — cos ky) the d-wave gap function and i?k = 
\/£,i^ + A^ the Bogoliubov quasiparticle energy. To facil- 
itate the derivation, we define two correlation functions 
L and M : 



Lu{p,z)^ ^ Gii{k + sp,z)A{k,uj), 



s=±,k 



M„(p,z)= ^ Gi2(k-|-sp,z)B(k,cj), 



(3) 



s=±,k 



where A(k, w) = — ImGii(k, aj"'")/7r is the single-particle 
spectral function and i?(k,aj) = — ImG2i(k, a;"'")/7r. It is 
straightforward to find that the ESS can be expressed as 

S{p,u) = Re[L^ip,Lo+) - AUp,u+)]. (4) 

On the other hand, the JDOS is given by 

J(p,c^)=^^(k + p,c^)A(k,w), (5) 
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and can simply be expressed in terms of : 

J{P,^) = --;^'^niL^{p,u}+). (6) 
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By examining Eqs. and (O, one can readily find 
that ^(pja;) and J(p,a;) are respectively related to the 
real and imaginary parts of the correlation function 
L^(p,aj+). Therefore they show singular behaviors on 
the common curves/surfaces in general, which stem from 
the interplay of the poles of the Green's function and 
spectral function in Eq. ([S]). This argument is verified in 
the following detailed discussions. 

To look into the ESS and JDOS of d-wave supercon- 
ductors more specifically, we now attempt to derive the 
explicit formulae of L^(p, z) and Mu:{p,z). The spec- 
tral functions A and B in Eq. ([3]) restrict that the mo- 
menta k's, which contribute to Li^ and AI^j, are de- 
termined by the four banana-shaped CCE's given by 
E]i ~ \lu\. Taking into account the d-wave structure 
of the gap, we here use the Dirac cone approximation 
around the four gap nodes ni^2,3,4 = (i??, where 
the index of each node is the same as the quadrant it 
locates. For instance, near the first node rii = (?7,?7), 
Ak sa UAfci, Ck ~ vph and sa v^(wAfci)^ + (wFfc2)^, 
where fci^2 = (k — ni) • (e^: =F ej,)/\/2 in the local coordi- 
nate system of the first node. The fci (^2) axis is along 
the direction of increasing Ak (^k)- If we perform the 
integration over k in Eq. (jS]) around the four gap nodes, 
analytical results can be obtained for the three cases: 1) 
p is small such that k and k -I- p sit at the vicinity of the 
same node, with q-j depicted in Fig. [TJa) being belong 
to this case; 2) p — (n^ — rii) is small such that k and 
k -|- p sit around the first and third nodes respectively, 
with q3_4 being belong to this case; 3) p — (n2 — ni) is 
small such that k and k -I- p sit around the first and sec- 
ond nodes respectively, with qi,2,5,6 being belong to this 
case. 

For the first case, after tedious derivations, we obtain 
L^(p,^+) = -;^-i— /+(p,L.+), (7) 
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/J(P,^+), 



, L,^-{Ep/2f±c,^cos^e , 1 
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where the angular factor cos 9 = v^pi/Ep. Substituting 
Eqs. © and (g]) into Eqs. ^ and ©, we have 

( 1 uj^ cos^ 8 01 ]\ ^ 

S{p,uj) = <i ^ '^"Ayp .e;p^e2_4^2 ' I I P' (^Q) 



0, 2\u\ > Ep, 



and 



J(p,a;) 



0, 2\uj\ < Ep, 

1 L^^-(gp/2)^+t^^ cos'' 



2\u\ > Ep. 



It can be seen clearly that both the ESS and JDOS 
show inverse-square-root divergence when approaching 
the same characteristic curve governed by Ep = 2\uj\, 
which will be termed as envelope in the following. 
The singular parts of both functions originate from the 
same complex function as shown in Eq. The ESS 
is nonzero (zero) outside (inside) the envelope while the 
JDOS is just the opposite. The vertices of this elliptic 
envelope are just the octet vectors iqy corresponding 
to 9 — 0,7r, on which both the ESS and JDOS reach 
the largest amplitude of divergence. This angular depen- 
dence of the ESS stems from the coherent factors of the 
d-wave superconductors as found in Ref. @ and here we 
find that the JDOS has the same angular dependence. 

For the second case, p' = p— (na— ni) is small. We find 
that the ESS and JDOS have the same analytic expres- 
sions as Eqs. pO| and (fTTj) if the corresponding numera- 
tors are respectively replaced by uj sin 9'{uj sin 9' — Ep/ /2) 
and (wsin^' — Ep'/2)'^, where sin6'' = VFp'2/Ep'. There- 
fore both quantities diverge on the envelope Epf = 2\uj\. 
However the angular behaviors of the ESS and JDOS 
are distinct from each other and from their behaviors 
in the first case. Accordingly both the ESS and JDOS 
do not reach their largest values on the vertices of the 
elliptic envelope i.e. ±q4 {9' = 0,7r). Another octet vec- 
tor qs locates on the envelope center where the JDOS 
has (5-function divergence, in agreement with the ARPES 
datcfc^. However the ESS vanishes at qa , seemingly incon- 
sistent with the FT-STS observation. The reason lies in 
the fact that the ESS reaches its strongest divergence at 
a vector of the envelope with 9' = — sgn(a;)7r/2, which is 
closest to q3 and therefore corresponds to the observed 
QPI peak as expected. 

For the third case with p near n2 — ni, we find again 
that both the ESS and JDOS diverge on certain common 
curves, although the explicit formulas for both quantities 
are rather involved than the above two cases. In order to 
capture the main physics without lengthy derivations, we 
propose an envelope scenario in the following to first lo- 
cate these curves where both the quantities exhibit diver- 
gent behaviors and then obtain the relative amplitudes 
of the divergence. 

To introduce the envelope scenario, we examine the 
JDOS in detail, which is relatively easier to handle. The 
spectral function of pure d-wave superconductors can be 
expressed as 

A(k,ij) ^ uld{uj - Ek) + vlS{uj + Ei,), 



where the coherence factors are 



(12) 

= 5(1 + k) ^'^d 

— ^(1 ~ "^)- Here we only study the uj < case 
for comparison with the autocorrclatcd ARPES data-. 
Therefore, we have after combining Eqs. ([5]) and (IT2t for 
w < 0, 



f d^k 

Jip,uj)^ J vl^pvlSiuj + E^+p)5iuj + Ek)j^ 

Leaving aside the coherence factors, one can easily find 
that the k points in the first Brillouin zone (FEZ), which 



(13) 



(11) 



3 



give the largest contribution to the integral, are deter- 
mined from the two CCE governed by the two 6 func- 
tions. 
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FIG. 1. Illustration of the idea of the envelope model. Here 
for simplicity we choose t' — ^ = and Aq = 0.2t. (a) CCE's 
determined by i5k = with uj = —OAt. The normal state 
Fermi surface is also shown, (b), CCE-1 in the first quadrant 
of FBZ. Two k points (red and blue dots) are shown, (c), two 
curves determined by i?k+p = \^\ with red solid corresponds 
to k point 1 while blue dashed one to 2 in (b). (d), the 
envelope is shown in green solid line and central red dot. For 
schematic illustration we also draw in blue color a cluster of 
curves for several different k's. qr as shown in (a) is located 
on the tip of the elliptic envelope. 



isolated point, 

Ep = 2\uj\, Ep = 0. (14) 

Taking advantage of the envelope functions, we can eval- 
uate the JDOS's divergence behavior just on these curves 
by integrating Eq. (|13p . which is relatively simpler than 
calculating J(p, uj) for arbitrary p. We find that the rel- 
ative amplitude of the inverse-square-root divergence is 
proportional to cos^ 9 on the elliptic envelope and the 
divergence is a 5 function on the isolated point. As for 
the inter-CCE case between CCE-1 and CCE-3, the en- 
velope consists of an ellipse and an isolated point similar 
to the intra-CCE case except that the center is shifted 
by a vector — rii and the relative amplitude on the 
ellipse changes to (1 + sin 6')^. All these results are in a 
full agreement with former explicit formulas. 
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FIG. 2. Plot of the envelope function according to Eq. (|15p . as 
well as 4 octet momenta qi,2,5,6. (a), less anisotropic case 7 = 
v^/vf = 0.4 (Ao — OAt, u — —0.8t). (b) strong anisotropic 
case 7 = va/vf = 0.1 (Ao = O.lt, uu = -0.2t). Note that the 
anisotropy is determined by Ao/t. We choose different u just 
to have similar octet g's for this two cases. 



In Fig. [1] we show the CCEs given by E^ = |a;|, which 
consists of 4 banana-shaped closed curves and are de- 
noted as CCE-1, 2, 3, 4 hereafter. To better elucidate our 
idea, we first focus on the intra-CCE contribution (i.e. 
both k and k + p on the same banana) to the JDOS. 
We examine the CCE in the first quarter of the FBZ, 
i.e. CCE-1. Corresponding to each k on CCE-1, all p 
points satisfying -Bk+p = |w| form a closed curve denoted 
as J7k(p), which has exactly the same shape as CCE-1 
except that the center is shifted as shown in Fig. [He). 
Therefore, for all k's on CCE-1, ?7k(p)'s form a one- 
parameter family of curves whose envelope is plotted in 
Fig. md) shown as the green line and red dot. One can 
readily see that the envelope gives the locus on which the 
JDOS is most likely strong because it represents the po- 
sition where curves in the family touch or intersect most. 

Quantitative analysis can be obtained under the Dirac 
cone approximation. Accordingly CCE-1 is an ellipse 
whose semi-axis lengths are a = \u}\/va and b = \ijj\/v-p. 
The corresponding family of curves I7k(p)'s as shown 
in Fig. [^d) are determined from the following equa- 
tions, {pi,P2) = {acos9 — ki,bsin9 — with {ki,k2) = 
{a cos (pjb sin (p), where 9,(j) £ [0,27r). The envelope of 
the above family is composed of an ellipse and a trivial 



More importantly, the envelope analysis can be ap- 
plied to the inter-CCE case between CCE-1 and CCE-2, 
which is otherwise rather complicated. The family of 
curves J7k(p)'s satisfy, (^1,^2) = — (fepifcr) + {bcos9 — 
ki,asm9 — ^2), where fcp — \/2ti denotes the length of 
the Fermi vector at the nodes. The corresponding enve- 
lope is 



pi = —kp + a cos a 
P2 = — fcp + a sin a 



7^5^,(0) ±r^(a)J , 
s.y{a) ± 7^r^(a) 



(15) 



where r^(a) — (cos-^ a -I- 7^ sin^ a) and Sj{a) = 
(sin^ a -f 7^ cos^ a)^-'^/^ with 7 = b/a — va/vf is the 
anisotropic factor. This envelope is plotted in Fig. [21 
where a circle-like ["-|-" sign in Eq. (|15p ] and an astroid- 
like curve ["— " sign in Eq. P3)) ] centered at (— fcp, —kF)i^ 
is clearly seen. Also shown are the 4 wave vectors (qi,2,5,6 
as in Fig. [IJa)) according to the octet model. One of 
our main findings is that, the anisotropy of the quasi- 
particle excitation spectrum plays a significant role in 
bridging the JDOS picture and the octet model. We 
find that the octet vectors qi,2,5,6 approach the cusps 
(most divergent as shown later) of the astroid for suf- 
ficient large anisotropy of the Dirac cone, while depart 
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from the cusps of the astroid for the less anisotropic 
case, as shown in Fig. [5] This indicates that the effec- 
tiveness of the octet model^i"— is actually associated with 
the strong anisotropy of the quasiparticle excitations in 
d-wave cuprates. 

On the circle-like and astroid-like envelopes, we again 
find the inverse square root divergence of JDOS similar 
to the intra-CCE case, with the relative amplitude (the 
coefficient of the divergent terms) , 



FIG. 3. Comparison of analytical and numerical results of 
the JDOS in the first Brillouin zone. The parameters are 
t' = -0.3t, fi = -l.Ot, Ao = 0.2t, uj = -0.2t. Left: Plot of all 
the envelopes. The gray levels on the curves signal the relative 
amplitudes of divergence. Right: Numerical results according 
to Eq. (fT3)) with lattice size 400 x 400, where S{x) 

with Tjo = O.Olt. 
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with 



g± (a) = [1 ± jr^ (a) sin a] [1 + 75^ (a) cos a] 



(16) 
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According to Eq. ((T^ . the amplitude F^{a) on the as- 
troid diverges when approaching its 4 cusps at a = 
7r/4, 37r/4, 57r/4, 77r/4, coinciding with the octet vectors 
qi,2,5,6 ill the large anisotropic case as shown in Fig. [5] 
This observation agrees well with the autocorrelation of 
the ARPES data^. The effect of the coherence factors is 
embodied by g±{a) in Eq. (fT6|). inducing an additional 
angular dependence of the amplitudes. Accordingly the 
amplitudes on the astroid- and circle-like envelopes are 
highly uneven with one half is rather larger than the other 
half, breaking the four-fold symmetry possessed by the 
envelopes. Thus these two pronounced parts combine to- 
gether to form a duck-foot structure as shown in Fig. [31 
Note that in the same spirit, the relative amplitude of 
divergence of the ESS on these envelopes can also be 
obtained (not shown here as the quantitative numerical 
analyses of the ESS have been presented before^"—). 
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Combining the above analyses, we can plot all the en- 
velopes where J(p, cj) diverges in Fig. [Ha) with the rel- 
ative amplitude of divergence being taken into consid- 
eration. We find two cigars centered at the origin, four 
duck feet on the coordinate axes and four (semi)cigars on 
the zone corners as well as five red dots. The envelope 
analysis capture the essential structures and symmetries 
of the JDOS, in excellent agreement with the experimen- 
tal data^ as well as the JDOS by numerical integration 
of Eq. as shown in Fig. [2Ib). As for the peak po- 
sitions, q7 and q4 locate on vertices of the central and 
corner cigars respectively. qi,2,5,6 is on apexes of the 
duck feet, qa lies on the center of the corner cigars de- 
noted as red dot. All these features are in full consistent 
with the experimental observations^. 

In summary, the JDOS and ESS in rf-wave supercon- 
ductors have been investigated analytically, both being 
found to diverge on the curves of envelopes. Using the 
envelope scenario by us, we have shown clearly that the 
the octet vectors locate on the vertices of the envelopes 
validating the octet model for the c?-wave cuprates with 
strong anisotropy of the quasiparticle excitation spec- 
trum. The relative amplitude of the divergence of the 
JDOS has been obtained. Finally, we wish to pinpoint 
that the envelope scenario may also be used to investi- 
gate the JDOS of the iron-based superconductors, which 
may shed light on the Fermiology and pairing symmetry 
in this new material. 
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In the local coordinate system the center of the astroid- and 
circle-like envelopes locates at (/cfjAf). More generally, it 
is n2 — ni. 



